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ABSTRACT 

We present the first 3D prints of output from a supercomputer simulation of a complex as- 
trophysical system, the colliding stellar winds in the massive (>120 M Q ), highly eccentric 
(e ~ 0.9) binary star system ?/ Carinae. We demonstrate the methodology used to incorporate 
3D interactive figures into a PDF journal publication and the benefits of using 3D visualiza¬ 
tion and 3D printing as tools to analyze data from multidimensional numerical simulations. 
Using a consumer-grade 3D printer (MakerBot Replicator 2X), we successfully printed 3D 
smoothed particle hydrodynamics (SPH) simulations of t) Carinae’s inner (r ~ 110 au) wind- 
wind collision interface at multiple orbital phases. The 3D prints and visualizations reveal 
important, previously unknown ‘finger-like’ structures at orbital phases shortly after perias- 
tron ((f) ~ 1.045) that protrude radially outward from the spiral wind-wind collision region. 
We speculate that these fingers are related to instabilities (e.g. thin-shell, Rayleigh-Taylor) 
that arise at the interface between the radiatively-cooled layer of dense post-shock primary- 
star wind and the fast (3000 kms -1 ), adiabatic post-shock companion-star wind. The success 
of our work and easy identification of previously unrecognized physical features highlight the 
important role 3D printing and interactive graphics can play in the visualization and under¬ 
standing of complex 3D time-dependent numerical simulations of astrophysical phenomena. 

Key words: hydrodynamics - binaries: close - stars: individual: Eta Carinae - stars: mass- 
loss - stars: winds, outflows 


1 INTRODUCTION 


The supermassive binary star system Eta Carinae (77 Car) is fa¬ 
mous for the greatest non-terminal stellar explosion ever recorded 
( [Davidson & Humphreys|fl997[ l. In the 1840s, 77 Car became the 
second brightest non-solar-system object in the sky and ejected be¬ 
tween 10 and 40 M 0 , forming the dusty bipolar “Homunculus” 
nebula ( [Smith et al.|2003||Gomez et al .|20 lTT[ |StefTen et al.|2014| (. 
Amazingly, this did not destroy the star(s). Multi-epoch ground- 
and space-based observations obtained over the past two decades 
reveal that 77 Car is a colliding wind binary (CWB) with a current 


total mass >120 M 0 ^Hillier 

et al. 2001 2006) and a highly eccen- 

trie (e ~ 0.9), 5.54-yr orbit ( 

Damineli et al. 1997; Whitelock et al. 

2004; Damineli et al. 2008a 

b, Gull et al. 2009, Groh et al. 2010b 

Corcoran et al. 2010, Teodoro et al. 2012). 


Hillier et 


Because they are so luminous (Z/rotal <; 5x 10 6 L ( 

|al.|2001[ |2006| l, the stars in 77 Car have powerful radiation-driven 
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stellar wind mass outflows. The Luminous Blue Variable (LBV) 
primary component, 77 A, has one of the densest known stellar winds 
(M nA * 8.5x 10 -4 Moyr -1 , Voo a 420 kms '^Hillier et al.|20(jl 
|Groh et aL||2012a] ). The less luminous companion star, 77 B, has a 
much lower density, but faster, wind ( L»/ L 0 ^ 10 s -10 6 , ^ 

1.4x 10~ 5 M Q yr -1 , v M as 3000 kms -1 ; Pittard & Corcoran|2002 


X-ray emitting shocks 

Pittard & Corcoran 2002[ [Corcoran et al. 

2010; Hamaguchi et al. 

2014) and a wind-wind interaction region 


(WWIR) that is thought to be the source of numerous forms of time- 
variable emission and absorption observed across a wide range of 
wavelengths (see e.g. |Damineli et al.|2008b[ . 


Proper numerical modelling of 77 Car’s WWIR remains a chal¬ 
lenge, mainly because it requires a full three-dimensional (3D), 
time-dependent treatment since orbital motion, especially dur¬ 
ing periastron, greatly affects the geometry and dynamics of the 
WWIR. Three-dimensional hydrodynamical simulations of 77 Car 
show that the complex time-varying WWIR has a major impact 
on the observed X-ray emission (Okazaki et al. 2008; Parkin et al. 
[ 20 TT] [R usse 11|2013) , the optical and ultraviolet (UV) light curves 
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and spectra 

Madura [2010, Madura & Groh 2012;|Madura et al. 

2012, 2013, 

Clementel^^df2014 2015), and the interpretation 


of various line profiles and interferometric observables |Groh et 
|al.| 2010 a|b|[ 2012 a|b| >. While 3D hydrodynamical simulations have 
helped to increase substantially our understanding of the present- 
day 77 Car binary, we are limited by an inability to adequately visu¬ 
alize the full 3D time-dependent geometry of the WWIR. However, 
this is crucial if we are to thoroughly understand how and where 
various forms of observed emission and absorption originate. 

Most published figures of 3D simulations of 77 Car’s collid¬ 
ing winds (and CWBs in general) consist of 2D slices through the 
origin of the typically-Cartesian 3D simulation domain, with color 
showing a scalar quantity such as density or temperature (see e.g. 
Lemaster et al .|2007|[Qkazaki et al.|2008llParkin et al.|2009|lPittard| 


2009, Parkin et al.|2011| and |Madura et al.|2013| hereafter M13). 
Sometimes, 2D slices showing physical quantities in the two coor¬ 
dinate planes perpendicular to the orbital plane, in parallel planes 
above and below the orbital plane, and in planes at arbitrary an¬ 
gles relative to the orbital plane, are also provided (Lemaster et al. 
|2007| [Okazaki et al.|2008||Pittard|2009| and appendix B of M13 1 . 
While multi-panel figures showing 2D slices can be useful, the 
time-varying geometry of the WWIR (caused by orbital motion), 
combined with parameter studies of various stellar, wind, and or¬ 
bital parameters, can lead to large numbers of cumbersome 2D fig¬ 
ures. The amount of information that 2D slices can convey about 
an intrinsically 3D structure is also limited, which can make such 
2D slices difficult to interpret and understand. 

There have been attempts to provide 3D isovolume or isosur¬ 
face renderings of the WWIR from 3D hydrodynamical simulations 
of CWBs, for example, fig. 3 of |Lemaster et al.|(2007j i and figs. 2, 
13, and 14 of Pittard (2009). For the specific case of rj Car, a few 
3D isovolume renderings exist in e.g. Parkin et al,| (| 201 l} , |Madura| 
|& Groh|(2012| l, and |Groh et al,| ( |2012b| . Such 3D renderings are 
typically the exception rather than the rule though, and, as is the 
case with 2D slices, multi-panel figures are necessary in order to 
show the full 3D structure of the simulation results from different 
viewing perspectives and/or at different orbital phases. 

The predominance of 2D figures and animations in the liter¬ 
ature is clearly driven by the need to display 3D data in a classic 
paper-journal format. In this sense, the problem of 3D visualiza¬ 
tion of complex simulations and observational data is not limited to 
77 Car or CWBs. However, there is no real reason that researchers 
should be limited to 2D graphics when presenting their results in 
peer-reviewed publications. This is especially true since all major 
astrophysical journals are now published online. So-called ‘aug¬ 
mented articles’ ( [Vogt & Shingles 1 2013 { > are possible, in which 3D 
interactive models, images, sounds, and videos can be included di¬ 
rectly within an Adobe Portable Document Format (PDF) article. 
The inclusion of 3D interactive models in the astrophysics litera¬ 
ture, via methods such as those described in |Bames & Fluke| ( |2008| ), 
are slowly becoming popular, and several astrophysical journals 
now fully support the inclusion of 3D interactive figures. Two very 
recent examples of such 3D interactive figures are fig. 1 of Vogt 
|et al.| (|2014|), which presents a novel new way to classify galaxy 
emission lines via a 3D line ratio diagram, and fig. 5 of |Steffen| 
|et al.| ( ]2014| ), which presents a 3D interactive model of Eta Cari- 
nae’s bipolar ‘Homunculus’ nebula that was constructed based on 
detailed spectral mapping observations obtained with the ESO Very 
Large Telescope/X-Shootej^] 


While the use of 3D interactive graphics will likely prove to be 
extremely helpful to astronomers in their quest for understanding 
and discovery, one should always try to keep an eye on emerging 
technologies that may further aid the astrophysical research com¬ 
munity. One such technology that has been increasing in popularity 
across many different fields in recent years is additive manufactur¬ 
ing or ‘3D printing’. 3D printing has the potential to provide an 
entirely new method for researchers to visualize, understand, inter¬ 
pret, and communicate their science results. 

The use of 3D printing in the astronomical community is still 
very much in its infancy, but several prominent examples have ap¬ 
peared within the past year. The first is a program aimed at trans¬ 
forming Hubble Space Telescope ( HST ) images into tactile 3D 
models using special software and 3D printers, with the goal of 
communicating the wonders of astronomy to the blind and visually- 
impaired ( jChristian et al.||2014| . In the peer-reviewed literature 
there exists a 3D printable version of fig. 1 from |Vogt et al.|(20 1 4[ > 
(see their fig. 15) and a 3D printa ble version of th e Eta Carinae Ho¬ 
munculus nebula model by Steffen et al. 12014 Each of these is a 
unique illustration of how observational data can be used to develop 
a 3D printable model for increased understanding and communica¬ 
tion of complex concepts. However, to date there have been no pub¬ 
lished attempts to use 3D visualization and printing techniques to 
aid in the understanding and communication of complicated multi¬ 
dimensional output from detailed numerical simulations of astro- 
physical phenomena. 

In an effort to further demonstrate the benefits of using 3D 
visualization and 3D printing as tools to analyze output and com¬ 
municate results from numerical simulations, we present the first 
3D prints of output from a supercomputer simulation of a complex 
astrophysical system, the colliding stellar winds in the t] Car binary. 
Using a consumer-grade 3D printer, we print output from 3D time- 
dependent smoothed particle hydrodynamics (SPH) simulations of 
77 Car’s inner (r~ 110 au) wind-wind collision interface at multiple 
orbital phases. The main goal is to gain an improved understand¬ 
ing of the full 3D structure of 77 Car’s WWIR and how it changes 
with orbital phase. The 3D prints and visualizations reveal previ¬ 
ously unrecognized ‘finger-like’ structures at orbital phases shortly 
after periastron that protrude radially outward from the spiral wind- 
wind collision region. The success of our work helps highlight the 
important role 3D printing can play in the visualization and under¬ 
standing of complex 3D time-dependent simulations. 

In the following section < 0 - we describe our methodology and 
numerical approach, including the SPH simulations, the 3D visual¬ 
ization of the SPH output, and the generation of 3D printable files. 
Section [3] presents the results in the form of standard 2D images, 
pictures, and 3D interactive graphics. A brief discussion of the re¬ 
sults and their implications is in Section [4] Section [5] summarizes 
our conclusions and outlines the direction of future work. 


2 METHODS 

2.1 The 3D SPH simulations 

The hydrodynamical simulation snapshots we visualize correspond 
to specific phases (apastron, periastron, and 3 months after peri¬ 
astron) from the Case A (,M nA » 8.5x 10 ~ 4 Mgyr -1 ) and Case C 
(M nA 2.4 x 1(L 4 MoyU 1 ) small-domain (r = 10a = 155 au) 3D 


1 For other recent examples, see references in Vogt & Shingles (2013 


- See also http://www.nasa.gov/content/goddard/astronomers-bring-the- 
third-dimension-to-a-doomed-stars-outburst/ 
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3D printing p Car’s colliding winds 3 


SPH simulations of |M13| In all simulations, the primary star’s 
wind terminal speed is set to 420 km s -1 , while the companion 
star mass loss rate and wind terminal speed are set to M m = 
1.4 x 1CT 5 M 0 yr _1 and v ra = 3000 kms -1 , respectively. Addition¬ 
ally, the SPH particle mass is 5.913 x 10 24 g for the wind of 77 a and 
5.9 1 3 x 10 23 g (2.9565 x 10 24 g) for the wind of % in the Case A 
(Case C) r = 10a simulations. We refer the reader to |M13| and ref¬ 
erences therein for details on the SPH simulations and an extensive 
discussion of the results. 

The spherical computational domain size of radius r= 155 au 
is chosen in order to investigate at sufficiently high spatial resolu¬ 
tion the structure of 77 Car’s inner WWIRs, since the ‘current’ inter¬ 
action between the two winds occurs at spatial scales comparable to 
the semi-major axis length a ~ 15.4 au. The three orbital phases se¬ 
lected are representative of when the WWIR has its simplest (apas- 
tron) and most complex (periastron and 3 months after periastron) 
geometry. The snapshot at apastron, when the stellar separation is 
largest and orbital speeds are lowest, provides a reference for the 
nearly-axisymmetric conical shape of the WWIR during most of 
the orbital cycle. The periastron snapshot defines the distorted ge¬ 
ometry of the WWIR when the stellar separation is smallest and 
orbital speeds are greatest. The snapshot at 3 months after peri¬ 
astron (0 = 1.045) corresponds to a time when the WWIR has a 
distinct Archimedean-spiral-like shape in the orbital plane due to 
the rapid orbital motion of the stars around periastron {Okazaki et| 
|al.|2008HParkin et al,|201 l||Madura et al.|2012||2013| >. Simulations 
assuming two different 77 ^ mass loss rates are used to investigate 
how changes to the wind momentum ratio alter the WWIR opening 
angle, apex distance, 3D geometry, and dynamics. 

We use a standard xyz Cartesian coordinate system and set 
the orbit in the xy plane, with the origin at the system centre of 
mass and the major axis along the x-axis. The stars orbit counter¬ 
clockwise when viewed from along the +z-axis. By convention, pe¬ 
riastron is defined as t = 0 (0 = f/2024 = 0). Simulations are started 
at apastron and run for multiple consecutive orbits. 

2.2 Grid construction and density distribution 

SPH is a mesh-free method for solving the equations of fluid dy¬ 
namics {Monaghan|2005[ that is widely used in the astrophysical 
community. However, visualizing SPH data is far from straightfor¬ 
ward since the data are highly adaptive and unstructured, defined 
on a set of points that follow the motion of the fluid. Simple inter¬ 
polation to a uniform structured grid is often not an option since the 
grids are so immense in size they cannot be handled efficiently, or 
significant interpolation errors are introduced in areas of high parti¬ 
cle density (Linsen et al. 201 l) l. Due to these complications, several 
programs have been developed specifically for the visualization of 
SPH data. A popular, freely-available tool that allows for visualiza¬ 
tion of slices through the simulation volume, direct volume render¬ 
ing, and particle rendering is splash ( |Price||2007) ). Unfortunately, 
splash does not currently allow 3D isosurface extraction, which is 
required if we want to visualize 77 Car’s stellar winds and WWIR as 
solid 3D surfaces. A more recent program designed for the interac¬ 
tive visual analysis of SPH data is smoothviz (Linsen et al .|2011| >. 
smoothviz allows isosurface extraction and direct volume render¬ 
ing. but is currently limited to producing standard screenshots for 
use in research papers, as it uses OpenGL instead of graphics li¬ 
braries such as pgplot ( [Molchanov et al.|2013| l. 

Since our goal is to create 3D interactive figures and 3D print¬ 
able files, we employ an alternative approach. We choose to gener¬ 
ate a tetrahedral mesh from the SPH particle data in order to facil¬ 


itate easier visualization with standard software such as VisU0and 
ParaView]^] The generation of tetrahedral meshes from particle data 
has a long tradition, with the widely accepted results of Delaunay 
tetrahedrization dating to 1934 & Wang||2006[ [Linsen et al. 


2011 1 . We therefore employ the same methodology as Clementel 


et al. ( |2014[ [2015[ > to generate our unstructured 3D mesh. Using 
the SPH particles themselves as the generating nuclei, we tessel- 
late space according to the Voronoi recipe: all points in a grid cell 
are closer to the nucleus of that cell than to any other nucleus. The 
Voronoi nuclei are then connected by a Delaunay triangulation. 

We assign to the nucleus of each Voronoi cell the correspond¬ 
ing SPH quantities of particle mass, density, temperature, and ve¬ 
locity, computed using the standard SPH cubic spline kernel {Mom] 
agha n|T992] l. This helps ensure that each scalar variable visualized 
on our mesh closely matches that of the original SPH simulations, 
since the kernel samples a larger number of particles over a larger 
volume, resulting in quantities that are less affected by local dif¬ 
ferences in the SPH particle distribution (Clementel et al. 2015). 
Comparison with a direct visualization of the SPH density using 
splash ( |Price|2007) shows that this approach indeed matches well 
the density distribution of the original SPH simulations (see fig. 1 
of Clementel et al. |2015| . Fig. 2 of Clementel et al. ( [201 5^ shows an 
example of the unstructured mesh and number density at apastron 
for one of our 3D SPH simulations of 77 Car. 


2.3 Visualization 

When visualizing scalar variables on our unstructured mesh, we 
would ideally like to render quantities that are centred on the orig¬ 
inal Voronoi cells that compose our 3D grid. Unfortunately, the 
Voronoi cells consist of a series of irregular n-sided polygons, 
which makes their visualization quite complex. Instead, it is much 
more straightforward to visualize the corresponding Delaunay tri¬ 
angulation. In 3D, the Delaunay cells are tetrahedra, which can be 
visualized using standard visualization tools. Since the Delaunay 
cells are tetrahedra, the quantity we visualize is the average of the 
four vertices that define the tetrahedron cell (i.e. the average of 
the four Voronoi nuclei). This approach works well for visualiz¬ 
ing most physical quantities (e.g. temperature, density, velocity), 
and is suitable for our work. However, if neighbouring Voronoi nu¬ 
clei have values which are significantly different (i.e. by several 
orders of magnitude), this ‘volume-average’ approach may lead to 
tetrahedral-cell values that are difficult to understand and interpret 
(see |Clementel et al.|2015| for details). 

To help the reader better appreciate the 3D structure of 77 Car’s 
WWIR and the cavity carved within 77 a’s wind by 773 , we provide 
in Figs. [TJ§ three related visualizations (columns) of each SPH 
model. In each figure, the first column shows an arbitrary view 
looking down on the orbital plane, with the lower-density 77 B wind 
cavity opening toward (top row) or away from (bottom row) the ob¬ 
server. The top half of the 3D model (z > 0) is transparent in order 
to clearly show the orbital plane. This view provides a useful refer¬ 
ence for comparing the 3D results to 2D orbital-plane slices shown 
in earlier works such as lM13l 

Unfortunately, in a view such as that in the first column of 
Figs. [TJ[ 6 ] the fully rendered wind of t/b prevents one from seeing 
the complete 3D geometry of the cavity carved within r/^’s wind. It 
would thus be useful to visualize the modified wind of ija separate 


3 https://wci.llnl.gov/simulation/computer-codes/visit 

4 http://www.paraview.org/ 
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from the lower-density 773 wind. Doing this is straightforward since 
77 A and /7 b have very different M (Mr/J Mj; b * 60). Due to this M 
difference, our SPH simulations use different particle masses for 
each stellar wind. Therefore, using the SPH particle mass, we can 
isolate the pre- and post-shock 77 A winds and visualize them while 
keeping the entire pre- and post-shock r 7 b winds, and the top (z > 0) 
half of the model, transparent. Examples of this view are shown in 
the middle column of Figs. 

Since we are most interested in the 3D structure of 77 Car’s 
WWIR, we must find a way to also visualize it separately from the 
individual stellar winds. For the WWIR, we choose to visualize the 
thin, dense post-shock 77 A wind region that is located between the 
contact discontinuity (which separates the colliding wind shocks) 
and the pre-shock 77 A wind. In order to isolate this specific region, 
we use its unique density and temperature distribution. Radiative 
cooling of the post-shock 77 A gas increases substantially its density, 
by at least an order of magnitude (Parkin et al. 20TTJ[MT3}. There 
is thus a large difference in density between the pre- and post-shock 
77 A winds. Since the SPH simulations use constant, spherical mass 
loss rates, the density dependence with radius from the stellar sur¬ 
face in each pre-shock stellar wind in the simulations is roughly 
given by PsphericalW = M/[4m 2 v(r)], where v(r) = v M (l -R+/rf 
is the standard ‘beta-velocity law' (/3 = 1 for our simulations), with 
Voo the wind terminal velocity and R* the stellar radius l |M13j >. By 
defining a quantity 6 = PSPH/Pspherical- we can determine at each lo¬ 
cation the contrast in density between what is provided in the SPH 
simulations (psph)> and what the density at that location would be 
in an undisturbed spherical stellar wind. Within the pre-shock wind, 
8 ~ 1, while in the post-shock wind, 8 is greater than one, usually 
much greater. Therefore, we isolate the dense, post-shock 77 A wind 
(from here on referred to generally as the WWIR) by computing 
8 for 77 a’s wind and selecting only those regions with 8 > 2 and 
T = 10,000 K. Because the post-shock 77 A wind cools radiatively, it 
remains at the floor temperature T = 10,000 K set in the SPH simu¬ 
lations (M13fr . Considering only regions with T - 10,000 K ensures 
that we isolate the WWIR from the much hotter (T > 10 6 K) post¬ 
shock 77 B gas. The last column of Figs.[T|j6]illustrates a view identi¬ 
cal to that in the middle column of the figures, but with the addition 
of the 3D WWIR surface that exists above the orbital plane. 

2.4 Generation of 3D interactive figures 

So that readers can directly experiment with and see for themselves 
the full 3D structure of 77 Car’s WWIR, we augment this article by 
incorporating directly into the PDF 3D interactive figures of both 
our SPH simulation results and the final 3D print models. We em¬ 
bed our 3D graphics directly into the PDF document in order to 
facilitate easy direct sharing of the 3D results. Nearly every figure 
in this paper has a 3D interactive counterpart that can be accessed 
using the freely-available software Adobe Acrobat Reader v.S.(Q 
or above. Unfortunately, interactive 3D graphics in PDFs can cur¬ 
rently only be viewed using Adobe Reader. Other PDF viewers will 
display only the standard 2D images shown in each figure. 

The 3D interactive graphics allow the reader to fully rotate, 
zoom, and fly around each 3D model. This is a very efficient tool 
for revealing the structure of 77 Car’s WWIR. In some cases, the 
3D figures allow for the display (or not) of different components 
of the model (e.g. Figs. [8]|10| and 1 1 lf[l 3| >, providing the reader 
control over what he/she wants to see. Pre-defined ‘views’ to help 

5 http://get.adobe.com/reader 


guide the reader to specific orientations or features are also imple¬ 
mented, such as the orientation of the 77 Car binary on the sky as 
seen from Earth. We highly recommend after selecting and click¬ 
ing on a specific 3D interactive figure, that the reader right-mouse- 
click the model and select from the available drop-down menu the 
option “View in Floating Window”. This will open the interactive 
3D model in a small (although adjustable) side window that permits 
continued reading of the text and simultaneous viewing of other 
figures with minimal inconvenience. When in 3D interactive mode, 
right-click and select “disable content” to return to the original 2D 
figure. Numerous other options are available in the toolbar asso¬ 
ciated with each 3D figure and we encourage the reader to fully 
explore these. 

To create our interactive 3D graphics, and in an effort to en¬ 
courage others to use such 3D figures in their work, we rely on 
robust freely-available visualization software. The most difficult 
part is converting a particular 3D visualization into the U3D for¬ 
mat required for embedding within the PDF document ( {Barnes &| 
|Fluke|2008]|Vogt & Shingles|2013[ l. We start by using either VisIt 
or Para View to open and visualize our unstructured grid data. The 
choice of VisIt or ParaView is optional, and any suitable scientific 
visualization program may be used, provided it outputs the created 
data to the required format. Once we create our visualization, we 
export the model as either an OBJ or X3D file, depending on what 
we want to show. Both formats preserve the detailed geometry of 
the 3D models, but when converted to the U3D format, we find that 
X3D can preserve the color table used in a detailed scientific visu¬ 
alization, whereas OBJ sometimes does not. The downside to the 
X3D format is that the file size is usually larger than that of OBJ 
since the full detailed color information is being stored. The OBJ 
format is useful for situations where only the overall geometry, a 
solid surface color/transparency, and a small file size are needed. 

Once an OBJ or X3D file is created it can be read directly 
into the freely-available 3D mesh processing software MeshLa#’| 
where the 3D model can be adjusted and corrections applied if 
needed, and then exported directly to the required U3D format. An¬ 
other option is to read the OBJ or X3D file into the professional 
(but still free) 3D rendering and animation software BlendefQ In 
Blender, textures and colors can be improved, added, subtracted, 
etc., and additional objects or meshes can be inserted (or removed). 
Numerous visualization possibilities exist with Blender, and we 
find that it is more stable and easier to use than MeshLab. How¬ 
ever, Blender does not support exporting directly to the U3D for¬ 
mat. Nonetheless, one can create their final 3D model for visualiza¬ 
tion using Blender and export it as either OBJ or X3D. Then, one 
need only use MeshLab to quickly and easily convert the OBJ/X3D 
file to the U3D format. In order to create the best possible visual¬ 
izations, and in anticipation of creating 3D printable STL files (see 
Section |Z5| >, we employ Blender in this work, with a conversion 
to U3D using MeshLab. The above process for creating a U3D file 
may seem cumbersome, but we find that is in fact quite straight¬ 
forward. The interested reader that does not require free software 
can alternatively purchase programs such as PDF3DReportGei^] 
for the creation of USDs and 3D PDFs. 

Once a U3D file is available, incorporating it into a PDF can be 
done using either standard commercial software or, as astronomers 
generally prefer, the free typesetting package IATgX. The movie 15 


6 http://meshlab.sourceforge.net/ 

7 http://www.blender.org 

8 http://www.pdf3d.com/ 
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and media9 DTgX style files fully support the embedding of 3D 
annotations in PDF documents. Using and calling the movielS and 
media9 packages is incredibly straightforward and requires almost 
no more effort than inserting a standard 2D figure. For this paper, 
we use the movie 15 package. 

2.5 Generation of 3D print files and printing the results 

Having a 3D model and mesh, even one that looks nice on a com¬ 
puter screen, does not necessarily mean that it is suitable for 3D 
printing. Beside the requirement of converting the 3D model to the 
appropriate hie format for use with a 3D printer (generally the STL 
format), a 3D design that is to be 3D printed must meet a few basic 
requirements. These ensure that the 3D model prints correctly. 

A 3D design hie to be printed must be closed, or “watertight” 
as it is typically referred. All components in the 3D model should 
be connected to create a solid. Objects must typically also be man¬ 
ifold, i.e. have no edges that are shared between more than two 
faces. There should also be no parts of the model that have zero 
thickness. Individual and floating points, lines, and planes should 
be removed, since these do not have three separate directions nec¬ 
essary for printing. Surface normals also need to all be pointing in 
the same direction outward from the surface of the model. This en¬ 
sures that the 3D printer does not ‘confuse’ the internal and outer 
surfaces of the model. Finally, as may be obvious, the 3D model 
must have the appropriate physical dimensions to be printed on the 
specific printer of choice (i.e. your model must fit in the printer). 

While different 3D printers and 3D ‘slicing’ software may 
have more stringent requirements, we found in our work that the 
above were the minimal requirements necessary for a successful 
3D print of our models. Creation of a 3D printable STL file follows 
nearly the same procedure as that outlined above for generating a 
U3D file for PDF display. One key distinction though is that most 
3D printers are monochromatic, meaning they print using a single¬ 
color material per extruder. As such, the ‘color’ of the 3D model 
to be printed does not matter, since typically only one or two solid 
colors at most will be available. Color information is also (usually) 
not stored in the STL file to be printed, and the choice of color is 
more of a manual ‘real-world’ decision than something chosen in 
the file creation phase. Since we are mostly interested in the 3D ge¬ 
ometry and dynamics of the WWIR, the main concern is preserving 
the overall 3D geometry of the design to be printed. 

Thus, to create our 3D print files, we again visualize our SPH 
simulation data with VisIt or ParaView, but once satisfied with a 
particular visualization, we export our design in the OBJ file for¬ 
mat. In order to fit our 3D model in our printer, while at the same 
time preserving as much interesting detail as possible, we crop the 
outer spherical edge and choose to visualize and print only the re¬ 
gion extending to radius r = lax 108 au from the system centre- 
of-mass. This is a fairly small cropping of the model and, for the 
orbital phases of interest, no crucial information is lost. 

Once we create an OBJ file, we import it into Blender, which 
has useful tools for detecting and correcting non-manifoldness, 
mesh holes, loose objects, and inverted normals. We find that in 
most cases, the most prominent problem with 3D OBJ files ex¬ 
ported by VisIt or ParaView is that they have a large number of 
inverted normals. Luckily, these are easily corrected in Blender 
via a few simple mouse clicks. Once the normals are corrected, we 
use Blender to remove any floating points, lines, or faces that may 
lead to printing errors, and we close any open ‘holes’ that make the 
model non-watertight. Here, we must be careful to mention that a 
large ‘hole’ can exist in a model in a general physical sense and 


still allow that model to be 3D printable (see Results below). In¬ 
stead, what matters is that the entire surface of the 3D model itself 
be closed, with all edges/faces connected. 

To help guide the reader, using Blender, we add to our model 
two small spheres that are connected by a small thin cylinder, which 
represent the stars and illustrates their location and separation with 
respect to each other, the dense r]\ wind, and the WWIR. In or¬ 
der to make the individual stars more visible at the scale of our 
models, we have increased the radius of each sphere to be a fac¬ 
tor 3.5 times larger than the correctly-scaled stellar radius. Thus, 
the spheres depicting the stars in Figs. [S||T4| have radii equivalent 
to 210 and 105 R 0 for 77 a and rj b, respectively. The correct stellar 
separation to scale is used at each phase. 

We make one other final adjustment in Blender before export¬ 
ing our model to the STL format. This adjustment further ensures 
that our model fits the 3D printer and has a stable base on which 
to stand once it is printed. To provide an interesting scientific ref¬ 
erence point for anyone viewing our 3D printed models, we rotate 
them to the correct derived orientation that the r] Car binary has on 
the sky as seen from Earth, with an inclination i = 138°, argument 
of periapsis a> = 263°, and position angle on the sky of the orbital 
angular momentum axis of PA- = 317° (see Madura et al.| |2012| 
|M13| >. Once rotated, we remove a small portion of the bottom of 
the model so that it has a flat base. We are careful to not remove 
any of the WWIR or cavity carved within the 77 A wind by 77 B. Only 
a small portion of the undisturbed outer spherical r/\ wind is re¬ 
moved. This has no effect on our results or conclusions, but allows 
the 3D printed models to be placed on a flat surface and oriented to 
provide the viewer with an idea of how the system and WWIR ap¬ 
pear on the sky, assuming North is up and East is left (see Fig. m 

To print our 3D models, we use a consumer-grade MakerBot 
Replicator 2X Experimental 3D Printer, which has dual-extrusion 
capabilities. We import our STL hies into the freely-available Mak- 
erWar^ 3D printing software and create the X3G hies specihc to 
the MakerBot printer. Since our 3D models are incredibly com¬ 
plex and, in many cases, contain free-hanging unsupported edges, 
we print our models in one color using one of the MakerBot’s ex¬ 
truders, and use the second extruder to print dissolvable support 
material. Once printed, we place our model in a limonene bath, 
which safely dissolves away the support material. We use the high¬ 
est layer resolution possible when printing (100 microns), and a 
physical size that occupies nearly the entire available build volume 
(model diameter a 6 inches = 15 cm across its widest part). 

Each 3D printed model consists of two parts, joined by small 
metal pins. The bottom half of each model consists of the dense 
77A wind and the hollow r/s wind cavity, while the top half consists 
solely of the WWIR (dense post-shock 77 A wind region). The two 
pieces are separable to allow the viewer to see the internal regions 
of the cavity carved by r /b in 7 /a’s wind, or the WWIR by itself (see 
e.g. Figs. 1 10| and 1 13). Additionally, we add two small beads con¬ 
nected by a short pin to represent the individual stars and illustrate 
their orientation and separation. The radii of the beads and their 
separation is to scale with the rest of the printed model, although as 
described above, the stellar radii have been increased by a factor of 
3.5 to make the stars more easily visible. 

Finally, we note that each 3D printable STL file is attached 
to this article as supplementary online material. Thus, anyone with 
access to a suitable 3D printer can in principle print their own tj Car 
wind model(s). The STL format should be compatible with most, 


9 https://www.makerbot.com/desktop 
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if not all, 3D printers currently on the market. The STL file should 
be ‘ready to print' without additional modifications, but due to dif¬ 
ferences in printers and printing software, we cannot guarantee that 
absolutely no modifications are necessary. At the very least, the 
model will need to be adjusted to physically fit the 3D printer being 
used. We also have only tried printing the models on a 3D printer 
with dual extrusion and dissolvable support material, and are un¬ 
sure of the results of using a single-extrusion printer. We welcome 
readers with questions about 3D printing the models, the model cre¬ 
ation/design, and/or the file creation process to contact us directly. 


3 RESULTS 

3.1 The 3D visualizations 


Figs. |TJj6] present the results of our 3D visualizations of 77 Car’s 
inner WWIR. Embedded within each 2D figure is an interactive 
3D version of the model view shown in the last column. Figs, 
correspond to the higher M^ A Case A simulation, while Figs, 
correspond to the lower M nA Case C simulation. Case A is more 


1 

3 

4 

6 


representative of t/a’s current M (Ml3), while Case C is included 


mainly to investigate how changes to M nA affect the overall 3D 
geometry of the WWIR. Each model shown has a radius r = la » 
108 au, measured from the system centre-of-mass. 

Figures [T] and [4] show that at apastron, the WWIR and cav¬ 
ity carved by 773 in 77 a’s dense wind have the expected nearly- 
axisymmetric conical shape, with the WWIR opening angle in¬ 
creasing with decreasing M nA |mT 3|. However, one does not fully 
appreciate how dramatic the change in WWIR opening angle is un¬ 
til it is visualized in 3D. In Fig. |T] the WWIR is conical, whereas 
in Fig. [4] it is nearly a plane. Another interesting feature more 
clearly visible in the 3D renderings is the very slight distortion of 
the WWIR in the direction of orbital motion. This is due to the ad¬ 
ditional component to the wind velocities caused by the non-zero 
velocity of the stars about the system centre-of-mass. 

We also see in Figs. ^ and that the WWIR is not a clean, 
smooth surface. Rather, it is corrugated and contains many smaller- 
scale bumps and protrusions that arise as a result of various insta¬ 
bilities, such as non-linear thin shell and Kelvin-Helmholtz l |Parkin| 
|et al.|20iT] |M 13j . While somewhat difficult to see in the density- 
colored models of Figs. |T]and[4] Figs. [ 8 | and 1 1 l|clearly show that 
the WWIR in the Case C simulation appears to contain much more 
detailed structure, including large trenches that extend into and 
around the surface of the WWIR that faces the observer. We are 
unsure if this means that the Case C result is intrinsically more un¬ 
stable, or if the WWIR is simply better resolved as this SPH sim¬ 
ulation used more SPH particles (by a factor of ~ 1.5) and had a 
slightly better numerical resolution. 

Figures [2] and [5] illustrate the twisted 3D geometry of the 
WWIR at periastron. In both figures we find that the leading arm 
of the WWIR (near the apex) is highly twisted in the direction of 
orbital motion, whereas the spatially-extended outer remnant of the 
WWIR trailing arm, created before periastron, maintains some of 
its initial axisymmetric geometry. There is also more curvature of 
the WWIR apex in the Case A simulation than in the Case C sim¬ 
ulation. In Case C (fig. [5]>, the WWIR is still amazingly flat and 
planar in shape, although with a ‘twist’ of the apex and a spatial 
displacement of the leading edge of the WWIR further into the pre¬ 
shock 7 /a wind. 

One particularly interesting feature in the WWIR at periastron, 
visible in both cases, is the hole in the trailing arm near the 


WWIR apex. This WWIR hole is prominently located above (and 
below) the orbital plane, and is something not noticed in previous 
works that relied on 2D image slices through the orbital plane. This 
hole is a direct result of the fact that, at those locations where the 
hole exists, there is no longer a WWIR. We emphasize that it is not 
due to our choice of parameters used to isolate the WWIR. Plots of 
the temperature show that there is no shock-heated gas in regions 
where the hole exists. The temperature of material in the WWIR 
hole is sj 10 4 K. Plots of the 6 parameter also show no density en¬ 
hancement for material where the WWIR hole is located. Material 
within the WWIR hole has <5 < 1. 

Instead, because of the extremely high orbital eccentricity and 
embedding of 77 b within ^a’s dense wind at periastron, 77b’s out¬ 
wardly expanding wind is unable to collide with 77 a’s downstream 
wind. The wind of 77 B becomes trapped in specific directions at pe¬ 
riastron and there is no longer any wind-wind collision at certain 
locations in the trailing wind. As a result, the unshocked primary 
wind starts expanding and filling in the low-density cavity that was 
carved by 77b’s wind during the broad part of the orbit. The out¬ 
ermost WWIR-trailing-arm located downstream from the hole re¬ 
mains intact at periastron simply because of time-delay effects; the 
hole near the apex has not had enough time to expand and propa¬ 
gate downstream. 

The physical reasons for the formation of the WWIR hole are 
demonstrated in Fig. [7] which shows the density and wind speed 
in the orbital plane at periastron for the small-domain (r = 1.5a) 
Case A simulation of |M13[ with the wind velocity vectors (arrows) 
overlaid. The length of the arrows is proportional to the magnitude 
of the wind speed. If the reader clicks the figure, a short movie 
plays showing the evolution of the density and wind velocity in 
the orbital plane, starting at orbital phase 0.95 (~ 100 days before 
periastron), continuing through periastron and ending at phase 1.05 
(~ 100 days after periastron). 

The animation shows that leading up to periastron, there is a 
typical wind-wind collision. However, shortly before and at peri¬ 
astron, the stellar separation is small enough that, combined with 
radiative inhibition effects, 77 B can no longer effectively drive a 
wind toward 77 A l [M13fr . The dense wind of 77 A thus overwhelms 77 B, 
which can now only drive a stellar wind in directions away from 77 A. 
Fig .[7] shows that at periastron, all of the wind vectors point away 
from 77 A, and 77 B is deeply embedded within 77a’s extremely dense 
inner wind. Thus, there can no longer be a wind-wind collision in 
the trailing arm downstream from the stars (downstream in the di¬ 
rection opposite that of the orbital motion). 

The results of Fig. [7] demonstrate that the WWIR hole is not 
due to an inadequate numerical resolution of the simulations, since 
Fig.[7]uses the highest-resolution, smaller-domain simulation from 
|M13| This simulation uses 8 times the number of particles as the 
r = 10a simulations shown in the figures throughout the rest of this 
paper. The SPH particle mass is 5.9 1 3 x 10 23 g for the wind of 77 A 
and 1.971 x 10 23 g for the wind of 77 B in the r - 1.5a simulations. 
Although different SPH particle masses are used for the individual 
stellar winds since t/a’s mass loss rate is much higher than that of 
77 B, the r = 1.5a simulations are of adequate resolution to resolve 
both stellar winds, as demonstrated in |M13[ 

We find that the appearance and size of the WWIR hole does 
depend on the mass loss rate used for 77 A. The hole at periastron is 
larger and more prominent in Case A than Case C. As 77a’s mass 
loss rate is lowered, the hole forms later, closer to periastron, and 
is smaller in overall size after periastron. This is consistent with 
the change in wind momentum ratio and WWIR opening angle as 
Mrj A is lowered. The fast wind from 773 is able to better compete 
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Figure 1. 3D view of 77 Car’s WWIR at apastron for the Case A ( M nA = 8.5 x 10 -4 M G yr _1 ) SPH simulation. The top left panel shows an (arbitrary) view 
looking down on the orbital plane, with the lower-density 773 wind cavity opening toward the observer. The top half of the 3D model is transparent in order to 
clearly show the orbital plane. The top middle panel is the same as the first, but with the lower-density 773 wind made completely transparent in order to better 
show the 3D structure of the 773 wind cavity. The top right panel is identical to the top middle panel, but also includes the 3D surface of the WWIR that exists 
above the orbital plane. Panels in the bottom row are identical to those in the top row, but show a 180° rotated view. Color shows log density in cgs units in all 
panels. The locations of the stellar winds and WWIR are indicated. Click the image for a 3D interactive version of the model view shown in the last column 
(Adobe Reader® only. We suggest selecting “View in Floating Window” in the right-mouse-click drop-down menu). When in 3D interactive mode, right-click 
and select “disable content” to return to the original 2D figure. A small white sphere in the 3D interactive model marks the location of the companion star 773 . 
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Figure 2. Same as Fig. [7] but at periastron. Click image for a 3D interactive model (Adobe Reader® only). 
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Figure 3. Same as l ig.JT] but at - 3 months after periastron (tj> = 1.045). Click image for a 3D interactive model (Adobe Reader® only). 
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Figure 4. Same as Fag.jT] but for the Case C SPH simulation (M nA 


= 2.4 X 10 4 M e yr 1 ). Click image for a 3D interactive model (Adobe Reader® only). 



against 7ja’s wind in Case C, and thus able to better maintain the 
trailing arm of the WWIR at periastron. Lower M r)k also delay and 
alter any WWIR collapse that may occur around periastron (M13 1 . 
This is consistent with the simulation results of |Parkin et al.| ( j2011 1 , 
since their simulations assumed a primary star mass loss rate that 
is roughly half that which we use for the simulations in Figs. [TJ|3] 
[7] and[8||TT)| in which the WWIR hole is most prominent. 


Interestingly, there is no evidence of a WWIR hole around pe¬ 
riastron in the 3D t] Car simulations of |Parkin et al.| ( [2011 [ >. There 
are several reasons for this, discussed extensively in Section 3.1.1 
and Appendix A2 of |M13| (to which we refer the reader for details). 
The key point is that the stellar radius and wind-velocity law as¬ 
sumed for j/b, together with radiative cooling effects, have a big 
influence on determining whether there will be a cooling-transition 
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Figure S. Same as Fig.^] but at periastron. Click image for a 3D interactive model (Adobe Reader® only). 



Figure 6. Same as Fig.^] but at - 3 months after periastron (<f> = 1.045). Click image for a 3D interactive model (Adobe Reader® only). 


phase from adiabatic to radiative in 77 g’s post-shock wind around 
periastron. The slightly larger stellar radius and beta-wind velocity 
parameter used in our simulations lead to a significant reduction 
in the pre-shock t/g wind speed around periastron, which results 
in strong, rapid cooling of the post-shock 773 gas. Using the 773 pa¬ 
rameters of ]Parkin et al.| ( [20TT| , the reduction in 773 ’s wind speed by 
77 A is insufficient to cause 773 ’s post-shock wind to switch strongly 


to the radiative-cooling regime and cause a WWIR collapse. 773 ’s 
wind in the simulations of |Parkin et al.| ( |201 1) 1 is thus able to effec¬ 
tively maintain a wind-wind collision during periastron passage, 
even downstream, explaining why no hole is seen in the trailing 
arm of the WWIR in their simulations. Therefore, the existence of 
the WWIR hole around periastron depends strongly on the assumed 
stellar, wind, and orbital parameters of the 77 Car system. 
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Figure 7. Density (left) and wind speed (right) in the orbital plane at periastron for the small-domain (r = 1 .5a) Case A simulation of |M13| Wind velocity 
vectors (arrows) are overlaid on both plots. The length of the arrows is proportional to the magnitude of the wind speed. Click the figure to play a short movie 
showing the evolution of the density and wind velocity in the orbital plane. The movie starts at orbital phase 0.95 (~ 100 days before periastron) and ends at 
phase 1.05 (~ 100 days after periastron). The movie frame rate is set to 15 frames/second in order to better show the evolution of the wind velocity. 


Of the three phases studied, that at ~ 3 months after periastron 
(0 = 1.045, Figs. mm and|13fr contains the most fascinating 3D 
WWIR geometry. As expected based on previous 3D simulations, 
a spiral cavity is carved within the dense wind of 77 a by 77 b’s fast 
wind. However, the full 3D geometry of this cavity has never been 
visualized before, and we observe some interesting new features 
and phenomena that occur in both M, ;a simulations. 

First, as alluded to above, as the 3D geometry of the WWIR 
evolves with time following periastron, the hole created near the 
WWIR apex at periastron expands and propagates downstream 
along the WWIR’s old trailing arm. Thus, the new spiral WWIR 
created during periastron essentially has no trailing arm, and con¬ 
sists predominately of a leading arm and structures above and be¬ 
low the orbital plane. The hole created at periastron grows as the 
system moves through periastron and slowly eats its way into the 
remnant of the WWIR that was created during the broad part of the 
orbit. This effect is more prominent the higher the value of M^ A . 

Next, we see that the spiral wind cavity carved by 77 b within 
the back side of wind during periastron passage is much shal¬ 
lower than the cavity carved during the broad part of the orbit. This 
is due to the short amount of time 773 spends on the far side of 77 A 
during periastron. The spiral wind cavity and new WWIR formed 
during periastron have yet to expand outward to a size comparable 
to the remnant cavity on the apastron side of the system. 

The most surprising new set of features found at 0 = 1.045 
are the protrusions or ‘fingers' that extend radially from the spiral 
WWIR (see Figs. [3] [6] [TO] and|13[l. Detailed examination reveals 
that these fingers are actually tubes that consist of a thin shell of 
cold (T ~ 10,000 K), dense, compressed post-shock 77 A wind that 
is filled with hotter (T > 10 5 K) post-shock 77 B wind. The fingers 
penetrate into the unshocked 77 A wind expanding on the periastron 
side of the system and extend slightly above and below the orbital 
plane from the leading arm of the WWIR. They do not extend per¬ 


fectly vertically above or below the orbital plane (i.e. they are not 
at an angle of 90° with respect to the orbital plane). Instead, they 
point radially outward away from the stars. The fingers extend in 
the same direction that 773 ’s fast wind is able to collide with 77 a’s 
wind during periastron passage, and their location appears to be 
tied to the direction and speed of orbital motion around periastron. 

Analysis of the high-resolution r = 1.5rt simulation shows that 
the first fingers start to develop at orbital phase 0 ~ 0.993 (about 
2 weeks before periastron). Additional fingers continue to develop 
until 0 * 1.01 (about 3 weeks after periastron). Approximately two 
dozen noticeable fingers are present on the entire 3D surface of a 
spiral WWIR once periastron passage is complete (i.e. at 0 > 1.01). 
The total number and distribution of fingers varies from spiral to 
spiral due to the nature of the instabilities that form them. The 
spacing of the fingers is generally comparable to a few times the 
thickness of the dense WWIR shell, although in some locations the 
spacing is approximately equal to the thickness of the WWIR shell. 

The length of the fingers is also impressive, being larger than 
the stellar separation at this phase (i.e. > 7 au). Their diameter is on 
the order of several au. The fingers never completely vanish in our 
simulations, but slowly expand and cool, increasing in volume and 
moving outward until they eventually leave the computational do¬ 
main. The expanded fingers are noticeable in the outermost shells 
of compressed 77 A wind located to the left in the large-scale SPH 
simulations shown in figures 8. B5, and B7 of ]M13[ These outer ex¬ 
panded fingers are quite large, in some cases larger than the entire 
central binary orbit. The gas within them has also cooled adiabati- 
cally down to the floor temperature of the simulations (10,000 K). 

We also find that M nA has some interesting effects on the 3D 
geometry of the WWIR at 0 =1.045. The WWIR in Case C appears 
to have fewer intact fingers. The protrusions present in Case C also 
appear to be shorter compared to those in Case A. Curiously, while 
the Case C WWIR lacks protrusions, it has an abundance of holes. 
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We suspect that these holes in the WWIR are actually indicative of 
protrusions, ones in which the dense outer shell of post-shock 77 A 
wind that defines the surface of the protrusion/tube has a S < 2. This 
idea is supported by examination of the wind cavity carved within 
77 a’s dense wind (Fig. 1 1 2|>, which exhibits small well-defined tubu¬ 
lar cavities that extend into the outwardly expanding unshocked rj& 
wind. Moreover, the holes in the Case C WWIR all line up and 
point in the same direction as the protrusions that are present. 

The spiral WWIR and wind cavity in Case C are also notice¬ 
ably broader than in Case A (Fig.| 1 3[>. This is not surprising and is 
a result of the larger opening angle due to the lower M^ A in Case C. 
What is intriguing is that in Case C, the outermost part of the lead¬ 
ing arm of the ‘current’ WWIR between/around the stars is so large 
that it has caught up to and started to collide and merge with the 
remnant of the old trailing arm that formed long before periastron. 
In Case A at </> =1.045, the current WWIR has yet to reach the re¬ 
mains of the old trailing arm and there is still a large separation 
between it and the remnants of the trailing arm. 

3.2 3D printing results 

Figs. |8}fT3| present our 3D printing results, showing a direct com¬ 
parison between the actual 3D printed model and a 3D rendering of 
the model. In each 3D interactive model, the default starting view 
has the model oriented at the same position on the sky as the 77 Car 
binary ( [Madura et al.|2012fr , i.e. the starting view shows how an ob¬ 
server from Earth would see the system on the sky. Additionally, in 
the ‘Views’ menu of the 3D graphics toolbar are options to display 
only the modified wind of rj a, only the WWIR, or both together. 
For those unable to view the 3D interactive graphics, we include 
Fig. |T4] which shows how each model would appear on the sky to 
an observer on Earth. The STL files used for the prints are available 
as supplementary material in the online version of this article. 

As demonstrated by the figures, the 3D printed models repro¬ 
duce remarkably well all of the key features observed in the 3D in¬ 
teractive visualizations. We were pleasantly surprised that individ¬ 
ual features in the 3D unstructured meshes, such as protrusions and 
trenches, were faithfully reproduced. One can make out in some lo¬ 
cations the tetrahedral grid cells used in the mesh. Even the small 
protruding fingers that extend radially from the WWIR at 0 = 1.045 
were reproduced and remained mostly intact. However, we did have 
one delicate finger on the Case A model accidentally break off, 
which was simply just glued back on (Fig.|10|l. 

The ability to hold and inspect the 3D printed models provides 
a new perspective on the WWIR’s geometry and an improved sense 
of the scale of the different structures. One appreciates more just 
how large the WWIR is compared to the stars and stellar separa¬ 
tion. The 3D models are also useful for constraining the observer’s 
line-of-sight to the binary and help demonstrate why certain lines- 
of-sight are inconsistent with available observations. For example, 
as illustrated in Fig. |T4] at apastron, for the assumed orientation, 
our line-of-sight lies within the WWIR cavity, implying that any 
X-rays generated at the WWIR apex would be detectable to an 
observer at Earth. However, rotating the apastron model by ~ 90° 
or more places our line-of-sight through the dense, optically-thick 
primary wind, which would absorb any X-rays emitted from the 
WWIR apex. Thus, we may safely rule out such lines-of-sight. 

The above is just one very simple example of how the 3D 
print models can be used. More importantly, we find that the 3D 
print models are extremely useful as a visual aid to help explain 
to non-;/ Car experts, and even non-astronomers, the 3D geometry 
and dynamics of the binary and WWIR. The 3D prints are a useful 



Figure 8. Comparison between the 3D rendering (top row) and 3D printed 
model (bottom row) of the Case A SPH simulation at apastron. Columns 
present a (arbitrary) view of the WWIR with the observer looking into the 
r/B wind cavity (left), and a view of the opposite side with the >/b cav¬ 
ity opening away from the observer (right). Both views are looking down 
on the orbital plane. Click the image for a 3D interactive version of the 
model (Adobe Reader® only). Pre-programmed views are available under 
the '“Views” menu in the 3D model toolbar. These include the projection of 
the system on the sky as viewed from Earth (view "LOS”, North up, East 
left), the primary wind and r/s wind cavity only (view "PrimaryWind”), 
and the wind-wind collision region plus stars only (view “WWCR”). The 
physical diameter of the 3D printed models is approximately 6 inches 
(15.24 cm), as measured across the flat orbital plane through the centre 
of the model. At this scale, 1 inch (25.4 mm) corresponds to a distance 
« 35.5 au as 5.31 X 10 9 km. The locations of the stars, WWIR, and orbital 
plane are indicated in the 2D figure. When in 3D interactive mode, right- 
click and select “disable content” to return to the 2D figure. 



Figure 9. Same as Fig.|8] but at periastron. Click image for a 3D interactive 
model (Adobe Reader*only). 
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Figure 10. Left 4 panels'. Same as Fig. [8] but at a 3 months after periastron (tp = 1.045). Right 4 panels : Comparison between a 3D rendering of the Case A 
simulation (top row) and the 3D printed model (bottom row) for the two separable pieces that compose the model at tp = 1.045. The left column shows the 
bottom half of the model with the dense wind and t/g wind cavity. The right column shows the top half that consists solely of the WWIR. The locations of 
the stars, orbital plane, WWIR, WWIR hole, and WWIR fingers are indicated. Click the figure for a 3D interactive model (Adobe Reader® only). 



Figure 11. Same as Fig. [8] but for Case C. Click image for a 3D interactive 
model (Adobe Reader® only). 



Figure 12. Same as Fig. ITT] but at periastron. Click image for a 3D interac¬ 
tive model (Adobe Reader® only). 


tool for illustrating concepts, relationships, and properties that are 
not easily conveyed by 2D, and even 3D, graphics. 


4 DISCUSSION 

The new 3D features discovered in our results may have some in¬ 
teresting implications for observational diagnostics of tj Car and 
other highly-eccentric colliding wind binaries, such as WR 140. 
The hole near the apex of the WWIR that appears in the simulation 
snapshot at periastron will affect the generation of the shock-heated 
gas responsible for t] Car’s observed time-variable X-ray emission 
l |Corcoran et al.|2010l|Hamaguchi et al.|2007||2014| >. Obviously, in 
areas where there is no wind-wind collision, there can be no shock- 
heated gas, and thus no thermal X-ray emission. Our 3D models 


imply that at periastron, there should be no thermal X-ray emis¬ 
sion from along the trailing arm. Furthermore, as M j;a is lowered, 
it should take longer for the hot gas, and thus X-ray emission, to 
vanish from the WWIR trailing arm as periastron is approached. 

The WWIR hole at periastron also provides us with informa¬ 
tion about the amount and type of material in line-of-sight at that 
time. Since there is no longer a WWIR directly in line-of-sight at 
periastron, the column density of material between us and the stars 
is dominated by unshocked primary wind material that is flowing 
to fill the wind cavity carved by t ]b during the broad part of the 
orbit. Interestingly, this situation makes ‘wind-eclipses’ by 77 A of 
various observed features (e.g. X-rays, Corcoran et al.|2010| and 
Hen emission, Teodoro et al. 20121 easier to achieve at periastron 
for two reasons. First, since there is no longer a trailing arm to the 
WWIR, the overall size of the WWIR and volume of shock-heated 
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Figure 13. Same as Fig.llOl but for Case C. Click the figure for a 3D interactive model (Adobe Reader® only). 


gas at periastron is much smaller, and therefore easier to eclipse. 
Second, the lack of a trailing arm to the WWIR allows the dense 
77 A wind to expand and fill the 77 B cavity in line-of-sight, increasing 
the size of the eclipsing wind photosphere. 


The observational implications of the fingers that protrude 
from the WWIR at phases ~3 months after periastron are unclear. 
Unfortunately, if real, they are too small to spatially resolve, even 
with HST. Hot gas within the tubular fingers may produce X-rays, 
but the intensity of any such X-ray emission is likely to be small 
compared to that of X-rays generated near the WWIR apex. It 
is also unclear if the shocks responsible for producing the fin¬ 
gers could contribute to the Hen T4686 emission observed across 
77 Car’s periastron passage. Portions of the fingers located within 
77 A’s inner He + zone may be able to produce a small amount of Hen 
emission if the shocks produce the required He + -ionizing photons. 
However, this is very speculative without more detailed modeling. 


It is difficult to constrain the observational implications of the 
WWIR fingers at this time because we lack a thorough understand¬ 
ing of their physical origins and properties. We currently speculate 
that the fingers arise as a result of instabilities at the interface be¬ 
tween the two colliding wind shocks, which undergo rapid complex 
changes around periastron due to the high orbital eccentricity and 
changing wind directions. We point out that during the creation of 
the spiral WWIR, a typical wind-wind collision does not take place, 
due to the much faster orbital motion of the stars. The leading arm 
of the spiral WWIR is formed not by a head-on collision of two 
spherical winds, but rather by the fast receding 77 B wind colliding 
with the much denser and slower receding 77 A wind. This collision 
still produces a pair of shocks, but the situation is now more analo¬ 
gous to that which occurs when a fast stellar wind interacts with a 
surrounding slower moving circumstellar shell (e.g. |Toala & Arthur] 
|2011( |van Marie & Keppens|201~2} |, with the post-shock 77 A wind 
forming a thin, dense shell via radiative cooling, and the post-shock 
773 wind remaining hot and cooling adiabatically (see Fig.[7j. 


Because of the high density contrast between the stellar winds, 
and because the faster 773 wind is pushing into and accelerating the 


slower 77 A wind, different instabilities, including thin-shell ( |Vish- 
niac| 1983] 1994 and Rayleigh-Taylor (RT|Young et al.|2001||ToaIa 


& Arthur|2011| van Marie & Keppens|20 12*^ are expected to arise 

and greatly distort/disrupt the spiral WWIR. Features very similar 


to the fingers that we observe in our models are visible in the wind- 
wind hydrodynamic simulations of Toala & Arthur (20 TT| (see e.g. 
their fig. 12). The protrusions that arise in the simulations of |Toala| 
|& Arthur] ( |2011 1 ) are a result of the thin-shell and RT instabilities. 
While not quite of the same magnitude as what we observe in our 
simulations, the 3D binary colliding wind simulations of |Pittard| 
< |2009[ their fig. 10) and |Parkin et al.|p011| their figs. 12 and 13) 
show somewhat similar features shortly after periastron, wherein 
the lower-density secondary wind penetrates into the denser pri¬ 
mary wind due to thin-shell instabilities. Since |Pittard| ( |2009} and 
Parkin et al. (2011) present mostly 2D slices of their simulations, 
we cannot be completely sure that the protrusions we observe are 
the same as the phenomenon shown in their figures. However, one 
reason such features may be stronger in our simulations is because 
our stellar winds have a larger density contrast than those in|Pittard| 
l [2009) , |Parkin et al.| ( [ 20 TTT ). 


We find no obvious dependence of the fingers’ properties on 
the resolution of the simulations in |M13| with the exception that 
higher-resolution simulations appear to better resolve the instabil¬ 
ities in the WWIR, possibly leading to more fingers. Our interpre¬ 
tation that the fingers are due to strong instabilities that form in the 
WWIR around periastron makes qualitative sense, but we note that 
standard SPH schemes are notorious for under-resolving shocks 
and instabilities jAgertz et al.|2007||Price|2008) . Thus, our results 
should be interpreted with caution until a more detailed analysis 
is performed. Simulations using a grid-based method will help to 
determine if the fingers are an artifact of the SPH scheme. Such 
simulations will be the subject of future work. However, the fact 
that the hydrodynamic simulations of |Toala & Arthur \2Q 11 |,[van| 
|Marle & Keppens| ( |2012) , |Pittard| ( |2009^ , and |Parkin et al.|(2011 1 
were performed using grid-based methods and produced qualita¬ 
tively similar phenomena supports the idea that the protrusions we 
observe are real. The timing and location of the fingers’ appearance 
also makes physical sense on the grounds that they appear only dur¬ 
ing the rapid periastron passage, and occur only in the directions in 
which the fast wind of r/s strongly collides with the slow, dense 
receding wind of r/\. 
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Figure 14. Comprehensive comparison between the 3D renderings and 3D printed models, with all models oriented as the system would appear on the sky to 
an observer on Earth (North is up. East is left). Rows (from top to bottom) present the three different orbital phases. The first two columns show, respectively, 
the 3D rendering and the 3D printed model, for Case A. The last two columns are the same as the first, but for Case C. 


5 SUMMARY AND CONCLUSIONS 

In this work we present the first 3D prints of output from 3D SPH 
simulations of a dynamic astrophysical system. We demonstrate the 
methodology used to incorporate 3D interactive figures into a PDF 
document and the benefits of using 3D visualization and 3D print¬ 
ing as tools to analyze data from multidimensional numerical simu¬ 
lations. In this paper, we investigate the 3D structure of the WWIR 
in the the innermost regions of the 77 Car binary system. Below we 
summarize our most important results. 

1. Several features in the overall geometry of the WWIR (e.g. 
changes in the opening angle due to different M nA , distortions 
of the cone cavity due to the orbital motion, and irregularities 
on the WWIR surface due to various instabilities) are fully and 
more easily appreciated with the help of 3D visualization. 

2. The inability of 773 ’s wind to collide with 77 A’s downstream wind 
produces a large hole in the trailing arm near the WWIR apex 
at periastron. The size of this hole is directly connected with the 
wind momentum ratio, and therefore the WWIR opening angle. 
As expected, the higher M n/K in Case A causes a larger hole than 
that in Case C. After periastron, the hole expands and propagates 
downstream along the vanishing WWIR trailing arm. 

3. The faster orbital motion during periastron passage produces a 
much shallower spiral wind cavity on the back side of 77 a’s wind 
compared to the nearly-axisymmetric conical cavity carved dur¬ 
ing the broad part of the orbit. 

4. The 3D models present new ‘finger-like’ features at phase 1.045. 
These protrusions extend radially, above and below the or¬ 
bital plane, outward from the spiral WWIR for several au. The 


WWIR in Case C exhibits large holes together with the pro¬ 
trusions. The presence of tubular cavities carved in r]\’s dense 
wind support the idea that these holes are actually indicative of 
protrusions with S <2. 

5. We speculate that the newly-identified finger-like protrusions 
are a result of thin-shell, RT, and other instabilities that arise 
where the receding fast 77 B wind collides with the dense, reced¬ 
ing 77 A wind. Future simulations using grid-based methods are 
needed to confirm the existence of the fingers and determine 
their physical origin and properties. 

6. The Case C model shows that at phase 1.045, the outermost part 
of the leading arm of the ‘current’ WWIR reaches and starts to 
merge with the remnant of the old trailing arm formed before 
periastron. In contrast, for Case A, there is still a large gap be¬ 
tween the ‘current’ WWIR and the remnants of the trailing arm. 

7. 3D models can be used to better visualize and constrain, for 
example, the observer's line-of-sight to the binary system, an 
important parameter to correctly interpret and model available 
observations. The 3D prints are also very useful for conveying 
complex ideas to non-experts. 

We demonstrate in this paper how software that generates 3D 
models, together with interactive 3D graphics in PDFs, can be 
used to produce publication-quality, scientifically instructive fig¬ 
ures. Moreover, we show that 3D printed models reproduce ex¬ 
tremely well the key features observed in the 3D interactive vi¬ 
sualizations. Even if we have only touched on the possible appli¬ 
cations of 3D printed models as a tool, this work helps highlight 
the important role 3D printing can play in understanding complex 
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time-varying astrophysical systems. The ability to physically inter¬ 
act with the 3D models provides a completely new way to visu¬ 
alize, analyze, understand, and disseminate such 3D simulations. 
3D print models are also extremely useful to show and explain to a 
non-expert or non-scientist the 3D geometry and dynamics of nu¬ 
merical simulations of astrophysical phenomena. Thus, 3D printing 
and visualization have the potential to improve the astrophysical 
community's ability to convey advances in our disciplines to the 
wider public, providing an opportunity for them to play a more ac¬ 
tive role in their learning by 3D printing various models. PDF is 
also the most widely-used, self-contained electronic document for¬ 
mat, implying that funding agencies, governments, and the public 
can easily interact with instructive, 3D representations of our work 
l |Bames & Fluke|2008f . We hope that this paper motivates others 
in the astrophysical community to pursue the use of 3D interactive 
visualization and 3D printing in their research and publications. 
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SUPPORTING INFORMATION 

Additional supporting information may be found in the online 
version of the article: 

3D Printable STL Files: 3D print files for the models shown in 
Figs.[§{T4] 

This paper has been typeset from a TgX/ IATgX file prepared by the 
author. 
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